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We consider impact parameter dependence of the polarized and unpolarized structure functions. 
Unitarity does not allow factorization of the structure functions over the Bjorken x and the impact 
parameter b variables. On the basis of the particular geometrical model approach we conclude that 
spin of constituent quark may have a significant orbital angular momentum component which can 
manifest itself through the peripherality of the spin dependent structure functions. 
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■ INTRODUCTION 

o 

The behaviour and dependence of the structure functions on the Bjorken x is among the most actively discussed 
subjects in the unpolarized and polarized deep-inelastic scattering. The particular role here belongs to the small x 
region where asymptotical properties of the strong interactions can be studied. The characteristic property of the 

1 ^ 1 low- a: region is an essential contribution of nonperturbative effects 0J^] and one of the possible ways to treat this region 

is the construction and application of models. Of course, the shortcomings of any model approach to the study of this 
• nonperturbative region is evident. However, one could hope to gain from these models an information which cannot 
be obtained by perturbative methods (cf. Q). Among possible extensions in these studies are the considerations of 
■ the geometrical features of the structure functions, i.e. the dependence of the structure functions on the transverse 
^ ' coordinates or the impact parameter. This subject is not new. The importance of the parton distributions in the 
r — . transverse plane was stressed in Q] and, e.g. a brief model discussion was recently given in |^]. This work is a revised 
and extended version of the latter one. As it has been demonstrated Q the 6-dependent parton distribution can be 
related to the Fourier transform of the off-forward matrix elements of parton correlation functions in the limiting case 
of zero skewedness. Impact parameter dependence would allow one to gain an information on the spatial distribution 
of the partons inside the parent hadron and the spin properties of the nonperturbative intrinsic hadron structure. The 
geometrical properties of structure functions play an important role under analysis of the lepton-nuclei deep-inelastic 
""^ ' scattering and in the hard production in the heavy-ion collisions. 

a: 

Qh' I- INTERPRETATION OF 6— DEPENDENT STRUCTURE FUNCTIONS AT SMALL x 

<D . 

We suppose that the deep-inelastic scattering is determined by the aligned-jet mechanism [EJ and consider the 
^-dependence of the structure functions along the lines used in The aligned-jet mechanism is an essentially 
nonperturbative and allows one to relate structure functions_with the discontinuities of the amplitudes of quark- 
5_j ■ hadron elastic scattering. These relations are the following 



q(x) = hm[F 1 (s,t)+F 3 (s,t)]\ t = (h Aq(x) = ^Im[F 3 (s, t) - F^s, t)]| t=0 , 5q(x) = hmF 2 (s, t)\ t=0 . (1) 

The functions Fi are helicity amplitudes for the elastic quark-hadron scattering in the standard notations for the 
nucleon-nucleon scattering. We consider high energy limit or the region of small x. 

The structure functions obtained according to the above formulas should be multiplied by the factor ~ l/Q 2 - 
probability that such aligned-jet configuration occurs [0. 

The amplitudes Fi(s.t) are the corresponding Fourier-Bessel transforms of the functions Fi(s,b). 

The relations Eqs. ([!]) will be used as a starting point under definition of the structure functions which depend on 
impact parameter. According to these relations it is natural to give the following operational definition: 

q(x,b) = hm[F 1 {x,b)+F 3 (x,b)], Aq(x, b) = ^Im[F 3 (x, b) - F x {x, b)], Sq(x,b) = hmF 2 (x,b), (2) 

and q(x), Aq(x) and 5q(x) are the integrals over b of the corresponding 6-dependent distributions, i.e. 

o 2 r°° o 2 r°° o 2 r°° 

qM = ^- bdbq(x,b), Aq(x) = -^- bdbAq(x,b), 5q(x) = / bdbSq(x,b). (3) 

7H2 Jq 7T X J 7T X J 
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The functions q(x, 6), Aq(x, b) and 5q(x, b) depend also on the variable Q 2 and have simple interpretations, e.g. the 
function q(x, b, Q 2 ) represent probability to find in the hadron a quark q with fraction of its longitudinal momenta x 
at the transverse distance 

b±Ab, A6-1/Q 

from the hadron geometrical center. Interpretation of the spin distributions directly follows from their definitions: 
they are the differences of the probabilities to find quarks in the two spin states with longitudinal or transverse 
directions of the quark and hadron spins. 

It should be noted that the unitarity plays crucial role in the direct probabilistic interpretation of the function 
q{x 1 b). Indeed due to unitarity < q(x, b) < 1. The integral q(x) is a quark number density which is not limited by 
unity and can have arbitrary non-negative value. Thus, the given definition of the 6-dependent structure functions 
is self-consistent. Of course, spin distributions Aq(x,b) and Sq(x,b) are not positively defined. 

II. UNITARITY AND £>— DEPENDENCE OF STRUCTURE FUNCTIONS 

The unitarity can be fulfilled through the [/-matrix representation for the helicity amplitudes of elastic quark- 
hadron scattering. In the impact parameter representation the expressions for the helicity amplitudes are the following 

F li3 (a;,6) = U h3 (x, b)/[l - iU 1>3 (x, b)], F 2 (x,b) = U 2 (x, b)/[l - iU x {x, b)] 2 (4) 

Unitarity requires lmU\, 3 (x, b) > 0. The [/-matrix form of unitary representation contrary to the eikonal one does 
not generate itself essential singularity in the complex x plane at x — > and implementation of unitarity can be 
performed easily. 

The model which provides explicit form of helicity functions Ui(x,b) has been described elsewhere Q. A hadron 
consists of the constituent quarks aligned in the longitudinal direction and embedded into the nonperturbative vacuum 
(condensate). The constituent quark appears as a quasiparticle, i.e. as current valence quark surrounded by the cloud 
of quark-antiquark pairs of different flavors. We refer to effective QCD approach and use the NJL model ^ as a basis. 
The Lagrangian in addition to the four-fermion interaction £4 of the original NJL model includes the six-fermion 
U(1)a -breaking term £g oc K(uu)(dd)(ss) 0]. Transition to partonic picture is described by the introduction of a 
momentum cutoff A = A x ~ 1 GeV, which corresponds to the scale of chiral symmetry spontaneous breaking . 

This picture for a hadron structure implies that overlapping and interaction of peripheral condensates in hadron 
collision occurs at the first stage. In the overlapping region the condensates interact and as a result virtual massive 
quark pairs appear. Being released a part of hadron energy carried by the peripheral condensates goes to genera- 
tion of massive quarks. In another words nonlinear field couplings transform kinetic energy into internal energy of 
dressed quarks. Of course, number of such quarks fluctuates. The average number of quarks in the cosidered case is 
proportional to convolution of the condensate distributions D® ,H of the colliding constituent quark and hadron: 

N(s,b)~N(s)-Df®D*, (5) 

where the function N(s) is determined by a transformation thermodynamics of kinetic energy of interacting conden- 
states to the internal energy of massive quarks. To estimate the N(s) it is feasible to assume that it is proportional 
to the maximal possible energy dependence 

N(s)~ K (l-(x Q ))s/t/(m Q ), (6) 

where (xq) is the average fraction of energy carried by the constituent quarks, (ttt-q) is the mass scale of constituent 
quarks. In the model each of the constituent valence quarks located in the central part of the hadron is supposed to 
scatter in a quasi-independent way by the produced virtual quark pairs at given impact parameter and by the other 
valence quarks. When smeared over longitudinal momenta the scattering amplitude of constituent valence quark Q 
may be represented in the form 

(f Q (s,b)) = [N(s,b) + N-l]{V Q (b)), (7) 

where N = Nh + 1 is the total number of quarks in the system of the colliding constituent quark and hadron and 
(Vq(6)) is the smeared amplitude of single quark-quark scattering. In this approach the elastic scattering amplitude 
satisfies the unitarity since it is constructed as a solution of the following equation 

F = U + iUDF (8) 
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which is presented here in operator form. The function U(s, b) (generalized reaction matrix) — the basic dynamical 
quantity of this approach — is then chosen as a product of the averaged quark amplitudes 



u( S ,b)= H(f Q (s,b)) 



(9) 



)=1 



in accordance with assumed quasi-independent nature of valence quark scattering. The strong interaction radius of 
the constituent quark Q is determined by its Compton wavelength and the ^-dependence of the function (/q) related 
to the quark formfactor Fq(q) has a simple form (/q) oc exp(— TOq&/£). The helicity flip transition, i.e. Q + — ► Q_ 
occurs when the valence quark knocks out a quark with the opposite helicity and the same flavor ill . 

The explicit expressions for the helicity functions Ui(x, b) at small x have been obtained from the functions Ui(s, b) 
H by the substitute s ~ Q 2 /x and at small values of x they are the following: 



ninX 



Ui, a (x,b) = U (x,b)ll + P h3 (Q 2 )m Q Vx~/Q}, U 2 (x, b) = g 2 (Q 2 )-^- exp[-2(a - l)m Q b/Z]U (x, b) 



where 



Uq(x, b) — iUo(x 1 b) = i 



g{Q 2 )Q 
rriQy/x 



N 



exp[-Mb/Z]. 



(10) 



(11) 



a, a, 0, gf and £ are the model parameters, some of them in this particular case of quark-hadron scattering depend 
on the virtuality Q 2 . The meaning of these parameters is not crucial here; note only that ttlq is the average mass of 
constituent quarks in the quark-hadron system of N = Nh + 1 quarks and M is their total mass, i.e. M — X)<=i m *- 
We consider here for simplicity pure imaginary case. We need to keep the subleading terms in the expressions for 
Ui(x, b) and U${x, b) since the Aq(x, b) is determined by their difference. For U2(x, b) one can keep only leading term. 
Then using Eqs. (Q) we obtain at small x: 



q(x,b) 



goM) 
l + U (x,bY 



Aq{x,b) 



P-{Q 2 )m Q ^x- U {x,b) 



Q 



[l + U (x,b)Y 



(12) 



q 2 AQ 2 )m 2 n x 
5q{x, b) = yfW ^ 2 Q exp[~2(a - l)m Q 6/fl 



U (x,b) 



[l + U (x,b)}< 



(13) 



where /3_ (Q 2 ) — j3^(Q 2 ) — fii{Q 2 )- From the above expressions it follows that q(x, b) has a central ^-dependence, while 
Aq(x, b) and 5q(x, b) have peripheral profiles. Their qualitative dependence on the impact parameter b is depicted in 
Fig. 1. The function Aq(x,b) has a maximum located at 

b max (x) = — In '—]. 

M mQ^/x 



q(x,b) 



X 2 <X 1 




Aq(x,b) 




FIG. 1. b-dependence of the structure functions q(x,b) and Aq(x,b) at low-x 
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From Eqs.([l^[L3|) it follows that factorization of x and b dependencies is not allowed by unitarity and this provides 
constraints for the model parameterizations of structure functions which depend on x and b variables. Indeed, it is 
clear from Eqs.( ^2|JT^ ) that factorized form of the input "amplitude" Uo(x,b) cannot survive after unitarization due 
to the presence of the denominators. It is to be noted here that from the relation of impact parameter distributions 
with the off-forward parton distributions 0] it follows that the same conclusion on the absence of factorization is also 
valid for the off-forward parton distributions with zero skewedness. 

The following relation between the structure functions Aq(x, b) and Sq[x, b) can also be inferred from the above 
model-based formulas 

6q(x, b) = c(g 2 )^ exp(- 7 &)A g (*, b). (14) 

Thus, the function 5q(x, b) which describes transverse spin distribution is suppressed by the factors \Jx and exp(— 76), 
i.e. it has a more central profile. This suppression also reduces double-spin transverse asymmetries in the central 
region in the Drell-Yan production compared to the corresponding longitudinal asymmetries. 

The strange quark structure functions have also a more central 6-dependcncc than in the case of u and d quarks. 
The radius of the corresponding quark matter distribution follows from Eq. ( [l3|) and is the following 

R q (x)~ Yl XnQ2/x (15) 

while the ratio of the strange quark distributions to the light quark distributions radii is given by the corresponding 
constituent quark masses, i.e. for the nucleon this ratio would be 

R s ( x )/R q ( x )^(l + ^L)-\ (16) 

where Am = mg — vhq. 

Time reversal invariance of strong interactions allows one to write down relations similar to Eqs.(P) for the fragmen- 
tation functions also and obtain expressions for the fragmentation functions D^(z, b), AD^(z, b), 5D^(z, b) which have 
just the same dependence on the impact parameter b as the corresponding structure functions. The fragmentation 
function D^(z, b, Q 2 ) is the probability for fragmentation of quark q at transverse distance b ± A b (Ab ~ l/<5) into 
a hadron h which carry the fraction z of the quark momentum. In this case b is a transverse distance between quark 
q and the center of the hadron h. It is positively defined and due to unitarity obey to the inequality < D^(z, b) < 1. 
The physical interpretations of spin-dependent fragmentation functions ADj(i, b) and SD' q l (x, b) is similar to that of 
corresponding spin structure function. 



DISCUSSIONS AND CONCLUSION 



It is interesting to note that the spin structure functions have a peripheral dependence on the impact parameter 
contrary to central profile of the unpolarized structure function. It could be related to the orbital angular momentum 
of quarks inside the constituent quark. The important point what the origin of this orbital angular momenta is. It was 
proposed jL2j to use an analogy with an anisotropic extension of the theory of superconductivity which seems to match 
well with the picture for a constituent quark. The studies Q of that theory show that the presence of anisotropy 

leads to axial symmetry of pairing correlations around the anisotropy direction / and to the particle currents induced 
by the pairing correlations. In another words it means that a particle of the condensed fluid is surrounded by a cloud 

of correlated particles which rotate around it with the axis of rotation I, Calculation of the orbital momentum shows 
that it is proportional to the density of the correlated particles. Thus, it is clear that there is a direct analogy between 

this picture and that describing the constituent quark. An axis of anisotropy I can be associated with the polarization 
vector of current valence quark located at the origin of the constituent quark. The orbital angular momentum L lies 
along I. 

Spin of constituent quark, e.g. t/-quark, in the model used is given by the sum: 

JU = 1/2 = S Uv + S{ qq } + L{qq\ = 1/2 + S{ qq } + £{gg}. (17) 

In principle, Si qq \ and L^ qq y can include contribution of gluon angular momentum, however, since we consider effective 
Lagrangian approach where gluon degrees of freedom are overintegrated, we do not concern problems of the separation 
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and mixing of the quark angular momentum and gluon effects in QCD (cf . |Q ) . In the N JL-model |T(]] the six-quark 
fermion operator simulates the effect of gluon operator f^G^G^, where is the gluon field tensor in QCD. It is 
worth to note here that in general large gluon orbital angular momentum is expected to be almost canceled by gluon 
spin contribution p5[ . 

The value of the orbital momentum contribution into the spin of constituent quark can be estimated according to 
the relation between contributions of current quarks into a proton spin and corresponding contributions of current 
quarks into a spin of the constituent quarks and that of the constituent quarks into proton spin |r| : 

(AE), = (AC/ + AD)(A£)t/, (18) 

where (AE)p = S Uv + S{q q y. The value of (AE) p was measured in the deep-inelastic scattering. Thus, on the 
grounds of the experimental data for polarized DIS we arrive to conclusion that the significant part of the spin of 
constituent quark in the model should be associated with the orbital angular momentum of the current quarks inside 
the constituent one 

Then the peripherality of the spin structure functions can be correlated with the large contribution of the orbital 
angular momentum, i.e. with the quarks coherent rotation. Indeed, there is a compensation between the total spin of 
the quark-antiquark cloud and its orbital angular momenta, i.e. L^ q y = —S{q q y and therefore this correlation follows 
if the above compensation has a local nature and valid for a fixed impact parameter. 

The important role of orbital angular momentum was known long before EMC discovery p7| and reappeared after 
as one of the transparent explanations of the polarized deep-inelastic scattering data . Lattice QCD calculations 
in the quenched approximation also indicate significant quark orbital angular momentum contribution to spin of a 
nucleoli jl9| . It is interesting to find out possible experimental signatures of the peripheral geometrical profiles of the 
spin structure functions and the significant role of the orbital angular momentum. One of such indications could be 
an observation of the different spatial distributions of charge and magnetization at Jefferson Lab [^| . It would also 
be important to have a precise data for the strange formfactor. 
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